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IE 1207.4-5209: THE PUZZLING PULSAR AT THE CENTER OF THE PKS 1209-51/52 SUPERNOVA 
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ABSTRACT 

Second Chandra observation of IE 1207.4-5209, the central source of the supernova remnant PKS 1209-5 1/52, 
allowed us to confirm the previously detected period of 424 ms and, assuming a uniform spin-down, estimate the 
period derivative, P ~ (0.7-3) x 10~ 14 s s . The corresponding characteristic age of the pulsar, P/2P ~ 200- 
900 kyr, is much larger than the estimated age of the SNR, ~ 7 kyr. The values of the spin-down luminosity, 
E ~ (0.4-1.6) x 10 34 erg s" 1 , and conventional magnetic field, B ~ (2-4) x 10 12 G, are typical for a middle-aged 
radio pulsar, although no manifestations of pulsar activity have been observed. If IE 1207.4-5209 is indeed the 
neutron star formed in the same supernova explosion that created PKS 1209-51/52, such a discrepancy in ages 
could be explained either by a long initial period, close to its current value, or, less likely, by a very large braking 
index of the pulsar. Alternatively, the pulsar could be a foreground object unrelated to the supernova remnant, but 
the probability of such a coincidence is very low. 

Subject headings: pulsars: individual (IE 1207.4-5209) — stars: neutron — supernovae: individual (PKS 
1209-51/52) — X-rays: stars 



1. INTRODUCTION 

Radio-quiet compact central objects of supernova remnants 
(SNRs) have emerged recently as a separate class of X-ray 
sources (see, e.g., Pavlov et al. 2002, for a review). The na- 
ture of at least some of these sources, presumably neutron stars 
(NSs), remains enigmatic. They are characterized by soft, ap- 
parently thermal, X-ray spectra and a lack of observed pulsar 
activity and optical counterparts. One of the most important 
properties, which potentially allows one to elucidate the nature 
of these sources, is the periodicity of their radiation. Our previ- 
ous observation of IE 1207.4-5209, the central object of PKS 
1209-51/52 (=G296.5+10.0), with the Advanced CCD Imaging 
Spectrometer (ACIS) on board of the Chandra X-ray observa- 
tory allowed us to detect a period P = 0.424129 s (Zavlin et al. 
2000; Paper I hereafter), which proved that the source is indeed 
a neutron star. This source was discovered with the Einstein ob- 
servatory (Helfand & Becker 1984), 6' off the center of the 81' 
diameter SNR. From the analysis of radio and optical observa- 
tions of this SNR, Roger et al. (1988) estimated the SNR age 
to be 7 kyr, with an uncertainty of a factor of 3. The distance 
to the SNR is d = 2.1^;| kpc (Giacani et al. 2000). The X-ray 
spectrum of the central source can be described by a thermal 
model, with a blackbody temperature of 3-4- MK and a radius 
of 1-2 km, at d = 2 kpc (Mereghetti, Bignami, & Caraveo 1996; 
Vasisht et al. 1997; Zavlin, Pavlov, & Trumper 1998). Fitting 
the spectrum with a hydrogen atmosphere model, Zavlin et al. 
(1998) obtained an effective temperature of 1.4-1.9 MK and a 
radius of about 10-12 km. Mereghetti et al. (1996) reported up- 
per limits on the source flux in radio (< 0.1 mJy at 4.8 GHz), 
optical (V > 23.5), and 7-rays (< 1 x 10~ 7 photons cm" 2 s" 1 for 
£>100MeV). 

In this Letter we present the results of second Chandra ACIS 
observation, which allowed us to estimate the period derivative 
and other parameters of the pulsar. 



2. OBSERVATION AND DATA REDUCTION 

IE 1207.4-5209 was observed with Chandra on 2002 Jan- 
uary 5-6 with the spectroscopic array ACIS-S in the Continu- 
ous Clocking (CC) mode. This mode provides the highest time 



resolution of 2.85 ms available with ACIS at the expense of 
one dimension of image. The source was imaged on the back- 
illuminated chip S3, at a focal plane temperature of -120 C. The 
total duration of the observation was 31.6 ks. Time history of 
detected events does not reveal substantial background flares, 
so we do not exclude any time intervals from the analysis. 

The ID image of IE 1207.4-5209 in "sky pixels" is con- 
sistent with the ACIS Point Spread Function, and it is simi- 
lar to the image obtained in the previous observation, taken 
at almost the same roll angle of 65.2° (vs. 65.3° in the sec- 
ond observation). To obtain the source count rate, we extracted 
20,282 source-plus-background counts from a ID segment of 8 
pixel length centered at the source position (10 pixels contain 
20,588 counts). The background was taken from similar seg- 
ments adjacent to the ID source aperture. Subtracting the back- 
ground, we find a source countrate of 0.64 ±0.01 s" 1 , versus 
0.76 ±0.01 s -1 in the previous ACIS observation. The reduc- 
tion of the count rate can be attributed to a reduction of ACIS 
sensitivity at low energies. 

3. TIMING ANALYSIS 

The event times in the event file of a CC observation are the 
times when the event is read out at the chip node. To restore 
the actual event arrival times, we subtracted the times of charge 
transfer between the detection point and the readout node from 
the readout times and corrected the times for the satellite wob- 
bling (dither) and the Science Instrument Module motion using 
the approach described in Paper I. The corrected times were 
transformed to the solar system barycenter using the axBary 
tool of the CIAO package (v. 2.2). 

From a 3 -pixel segment centered at the source position, we 
extracted 18,237 counts in the 0. 3^4.0 keV energy range, where 
the background contribution is negligible (about 0.7%). We ran 
the Z ; 2 test (m is the number of harmonics included — see Buc- 
cheri et al. 1983) in the vicinity of the pulsation frequency de- 
termined in our previous observation. For m = 1, 2, 3, 4, and 
5, we found the Z 2 peak values of 62.2, 63.9, 73.0, 77.6, and 
79.3 at frequencies lower than a reference frequency 2.35777 
Hz by 5.2, 4.9, 2.4, 2.8, and 3.3 /jMz, respectively. The proba- 
bility to find such high peaks in a noise spectrum is extremely 
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Fig. 1. — Power spectra (upper panel) and frequency depen- 
dences of odds ratio (middle and bottom panels) around the 
pulsation frequency for the first and second observations (solid 
and dashed curves in the upper panel, respectively). Dotted 
lines show the median frequencies f° 2 found in the two obser- 
vations. 



low, ~ 10" 14 in a frequency range <; 10 ^Hz, which confirms 
our previous period detection. For a uniform comparison of the 
peak frequencies in first and second observations, we recalcu- 
lated the Z\ peaks for the first observation using exactly the 
same extraction segments, energy range, and time corrections 
as in the second one and obtained = 61.6, 61.8, 62.3, 62.7, 
and 65.6 (notice the smaller contributions from higher harmon- 
ics in first observation compared to second one, in which third 
and fourth harmonics are quite significant). The peak frequen- 
cies in the second observation are lower than in the first one: 
fx -f 2 =4.2, 2.9, 1.1, 1.9, and 2.2 ^Hz, for m = 1, 2, 3, 4, and 5, 
respectively. This result shows that the pulsar has slowed down 
in two years, although the scatter of the frequency shift is rather 
high. 

To evaluate the most plausible frequencies and the frequency 
shift, we employed the method of Gregory & Loredo (1996), 
based on the Bayesian formalism. This method uses the phase- 
averaged epoch-folding algorithm to calculate the frequency- 
dependent odds ratio O(f). This ratio specifies how the data fa- 
vor a periodic model of a given frequency / over the unpulsed 
model, and it allows one to find the corresponding probability 
distribution p(f) oc O(f) f' 1 (see Gregory & Loredo 1996, Pa- 
per I, for details). The profile O(f) calculated for the second 
observation appears to be much broader than for the first one 
(FWHM^ 10 yuHz vs. ~ 3 /uHz), with several narrow peaks 
of comparable heights (see Fig. 1). Exploring the O(f) de- 
pendence in various energy ranges, which show different rel- 




0.00 



ft-fi [MHz] 

Fig. 2.— Probability distribution for the frequency shift (in 
^Hz _I ). The vertical dash line shows the median shift. 



ative contributions of the harmonics of the pulsation frequency, 
we found that the separate peaks are associated with different 
harmonics (e.g., the highest peak in the lower panel of Fig. 1 
is associated with m = 1), in agreement with the different 
peak frequencies mentioned above. The median frequency of 
the probability distribution in second observation, for the 0.3- 
4.0 keV band, is f° = 2.3577625 Hz; it differs from the mean 
frequency by ^0.1 /iHz. The uncertainties of the frequency 
are (-3.2, +4.5), (-4.2, +8.0), and (-4.5, +10.1) ^Hz at 68%, 
90%, and 95% confidence level, respectively; the standard de- 
viation is ct/2 = 3.7 /zHz. We repeated a similar calculation 
for the data from the first observation (using the same extrac- 
tion segment and energy range as for the second one) and ob- 
tained f° = 2.3577699 Hz, with uncertainties of (-1.4, +1.2), 
(-2.5, +1.8), and (-2.9,+2.1)at68%, 90%, and 95% confidence 
level, respectively; <jf \ = 1.3 ^iHz. 

We calculated the probability density distribution for the fre- 
quency difference, A =f 2 -f u as p(A) = J p 2 (f) pi(f-A)df 
(see Fig. 2), where p\(f\) and p 2 (ji) are the probability den- 
sities obtained for the two observations. The median of this 
distribution, A° = -7.0 ^Hz, gives a shift somewhat smaller 
than the mean shift, / A p(A)dA ~ fi-f° = -7.4 ^Hz, but the 
difference is well within the uncertainties of A: (-3. 8, +4. 5), 
(-5.4.+7.1), and (-5.9,+8.3) /xHz at 68%, 90%, and 95% con- 
fidence level, respectively; cr\ = 3.9 /jHz. Because of the 
large widths of the p 2 (fi) and, consequently, p(A) distribu- 
tions, there is a non-negligible (albeit small) 5.6% probability 
that/2 >/i. 

Adopting A = -7.0^ | /xHz for the frequency shift and 
assuming a steady slowdown (lack of strong timing noise, 
glitches, etc) during two years between the observations, we 
obtain the following estimates for the time derivatives 



/ = -l.l^xl0 _ "Hz s" 



■■2.0 + _\jx lO-^ss" 1 



(1) 



The light curve extracted at / = f% in the energy range 0.3- 
4.0 keV (Fig. 3, right panel) reveals one broad pulse per period 
with a source pulsed fraction of f p = 8 ± 2%, similar to that 
obtained in first observation. Figure 3 indicates that in both ob- 
servations there is a considerable shift of the pulse phase, up 
to 0.4-0.6, between the energy bands 0.3-1.0 and 1.0-1.7 keV, 
perhaps accompanied by a change of the pulse shape. (Similar 
behavior has been observed in a number of middle-aged pulsars 
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— e.g., Ogelman 1995.) There is also some evidence that the 
pulse shapes are different in the two observations — e.g., in the 
0.3^1.0 keV band the pulse is closer to a sinusoid in first ob- 
servation while its shape becomes more "triangular" in second 
observation. 



4. DISCUSSION 

The detection of approximately the same period in the sec- 
ond observation confirms our previous detection and proves 
unambiguously that IE 1207.4-5209 is a neutron star. How- 
ever, the time derivative of the period appears to be surpris- 
ingly small — the corresponding characteristic age of the pul- 
sar, t c = P/(2P) = 340!?!$ kyr, exceeds the estimated SNR age, 
3-20 kyr, by a large factor, and even the most conservative 
upper limit, P < 5 x 10" 14 s s" 1 , yields an incredibly old age, 
t c > 130 kyr. Such values of P and P, as well as the corre- 
sponding values of rotation energy loss rate, E = 4ir 2 IPP~ 3 = 
1.0^0 7 x 10 34 /45 erg s" 1 , and conventional "magnetic field", 

B = (3Ic 3 PP/8tt 2 R 6 )^ 2 = 2.9!?;? x W 2 I% 2 Rf G, are typical 
for radio pulsars. 

This result could be easily explained if one assumes that 
IE 1207.4-5209 is merely a middle-aged field pulsar, unrelated 
to the SNR. If this pulsar were a background object at a dis- 
tance well beyond the SNR, it would be difficult to explain the 
very high (thermal) X-ray luminosity — e.g., L x « 3 x 10 34 erg 
s" 1 at d = 10 kpc cannot be emitted from radio pulsar polar caps 
because it exceeds E, and it cannot be emitted from the surface 
of a cooling NS because the observed temperature is too high 
for such an age. One, however, might assume that the pulsar 
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Fig. 3.— Light curves extracted in different energy ranges 
for the first and second observations (periods 0.4241296 s and 
0.42413093 s, respectively). 



is a foreground object, at a distance of several hundred par- 
sees, perhaps similar to the radio-quiet Geminga (r c ~ 300 kyr, 
d ~ 0.2 kpc) or (radio-bright) PSR B1055-52 (r c ~ 500 kyr, 
d ~ 1 kpc). In this case, the observed thermal-like pulsed X- 
ray radiation could be interpreted as emitted from hot polar caps 
with a size of a few hundred meters. The lack of detected ra- 
dio and 7-ray radiation might be due to unfavorable orientation 
of pulsar beams (albeit 7-ray beams are expected to be rather 
broad). To explain the lack of softer thermal radiation from 
the entire cooling NS surface (like that observed from Geminga 
and PSR B 1055-52), one could speculate that IE 1207.4-5209 
is older than the above-mentioned pulsars (e.g., r c ~ 1 Myr can- 
not be ruled out), so that its surface temperature is lower, and its 
thermal radiation is too soft to be observable in X-rays. How- 
ever, although the observed properties of IE 1207.4-5209 can- 
not exclude the hypothesis that it is a foreground pulsar, the 
probability to find a middle-aged pulsar so close to the line of 
sight towards the SNR center is very low. Therefore, we con- 
sider this interpretation unlikely. 

If IE 1207.4-5209 is indeed associated with the SNR, we can 
reconcile the pulsar's age with that of the SNR, assuming that 
the braking index n of the pulsar is large and/or its initial period 
P() was close to its current value. If the pulsar slowdown is 
described by the equation / = -Kf n , then its age, for a constant 
K,is 

(n-l)P V P 



1- ^ 



(2) 



To obtain r ~ 20 kyr (an upper limit on the age of PKS 1209- 
51/52), one can assume, for instance, n ~ 20-90 at Pq <C P, 
or Pq ~ 395—417 ms for n = 2.5. We are not aware of phys- 
ical models of pulsar deceleration which would give so large 
braking indices, although Johnston & Galloway (1999) argue 
that empirical braking indices of relatively old pulsars can be 
as large as 20-50. On the other hand, there is no physical limi- 
tations on the initial period — it can be very close to its current 
value if the constant K in the slowdown equation is small. 

It should be noted that IE 1207.4-5209 is not the only pul- 
sar that shows a characteristic age much older than the age 
of the SNR with which it is apparently associated. A sim- 
ilar example is the 7-s Anomalous X-ray Pulsar (AXP) IE 
2259+586 associated with the SNR CTB 109 — the SNR age, 
~ 3-21 kyr (Rho & Petre 1997; Parmar et al. 1998), is much 
younger than t c ~ 226 kyr (Kaspi, Chakrabarty, & Steinberger 
1999). Another example is the 65 ms rotation-powered X-ray 
pulsar AX J181 1.5-1926, whose period and period derivative 
were measured with ASCA (Torii et al. 1997, 1999), and as- 
sociation with the remnant Gl 1.2-0.3 of the historic super- 
nova A. D. 386 was strongly supported by Chandra observations 
(Kaspi et al. 2001a). The characteristic age of AX J1811.5- 
1926, about 24 kyr, is 15 times the true age. 

If the pulsar was born in the supernova explosion that created 
PKS 1209-51/52, its X-ray luminosity, L x w 1 x 10 33 erg s _1 in 
the 0.5-6.0 keV band (at d = 2 kpc), is 0.06-0.3 of the spin- 
down luminosity E. This means that IE 1207.4-5209 is differ- 
ent from AXPs, which show L x > E. On the other hand, the 
fraction of the spin-down luminosity detected in X-rays is too 
large to interpret the X-ray radiation as a nonthermal radiation 
or thermal polar-cap radiation emitted from a rotation-powered 
pulsar (in fact, the X-ray spectrum of IE 1207.4-5209 is incon- 
sistent with a nonthermal, power-law model — see Zavlin et 
al. 1998). Therefore, the original interpretation of the spectrum 
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as emitted from the surface of a cooling NS (Matsui, Long, & 
Tuohy 1988) remains most plausible, and the size of the emit- 
ting region and the effective temperature can be reconciled with 
those expected for a 10-kyrNS assuming that the NS is covered 
with a hydrogen or helium atmosphere (Zavlin et al. 1998). The 
magnetic field ~ 3 x 10 12 G, inferred from the estimated P as- 
suming a centered magnetic dipole and /?ns =10 km, may look 
too low to cause the temperature non-uniformity (Shibanov & 
Yakovlev 1996) required to explain pulsations detected from 
IE 1207.4-5209. However, we can speculate that the magnetic 
dipole is off-centered, which can give a much stronger local 
magnetic field at the same value of the magnetic moment. 

While estimating the period derivative fromjust two observa- 
tions, one should not forget that such an estimate implies a lack 
of glitches or other timing irregularities. Young and middle- 
aged radio pulsars show glitches with Sf up to 5 x 10~ 6 / 
(Lyne & Graham-Smith 1998). Even higher timing irregular- 
ities have been observed from some AXPs (e.g., Sf ~ 10 -4 / in 
IE 1048.1-5937 — Kaspi et al. 2001b). If, for instance,' our 
pulsar had a glitch Sf ~ 10~ 5 / just before the second observa- 
tion, the pre-glitch frequency would be lower than the measured 
one by ~ 20 /iHz, which would increase the frequency shift and 
lowered the NS age by a factor of 3. Finally, we cannot rule out 
that IE 1207.4-5209 is in a binary system with a long orbital 
period — our X-ray observations, each spanning about 8 hours, 
suggest that the period should be P or b ;> 40 hours, as no effect of 
the binary orbit is clearly observed in the data, and the optical 
observations were not deep enough to exclude the presence of a 
faint secondary companion. In this case, the X-ray observations 
give Doppler-shifted periods — e.g., a radial velocity of 10 km 
s" 1 corresponds to a shift of 79 /iHz, much larger than the ob- 
served one. Therefore, the above estimates of pulsar parameters 
should be taken with caution until more timing observations of 
the pulsar are carried out. 

To conclude, we measured the shift of the pulsar period be- 
tween two Chandra observations. The period and its apparent 
derivative do not contradict to the hypothesis that it is a young 
pulsar physically associated with the SNR PKS 1209-51/52. 
From the data available, we cannot completely rule out that 
that this is a foreground middle-aged pulsar, but this alterna- 
tive interpretation looks less plausible. The inferred properties 
of the pulsar warrant deep radio, optical, and 7-ray observa- 
tions of this object. To understand the nature of the pulsar, its 
timing behavior should be monitored in X-rays, which would 
allow one to find out whether it shows strong timing irregular- 
ities or it is in a binary system. Another clue to the nature of 
IE 1207.4-5209 would be provided by deep, high-resolution 
spectral observations of this object in X-rays, which could de- 
tect and identify spectral lines and evaluate the surface mag- 
netic field of the neutron star. 
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